Nonequilibrium conductance through a benzene molecule in the Kondo regime 
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Starting from exact eigenstates for a symmetric ring, we derive a low-energy effective generalized 
Anderson Hamiltonian which contains two spin doublets with opposite momenta and a singlet for the 
neutral molecule. For benzene, the singlet (doublets) represent the ground state of the neutral (singly 
charged) molecule. We calculate the non-equilibrium conductance through a benzene molecule, 
doped with one electron or a hole (i.e. in the Kondo regime), and connected to two conducting 
leads at different positions. We solve the problem using the Keldysh formalism and the non-crossing 
approximation (NCA). When the leads are connected in the para position (at 180 degrees), the model 
is equivalent to the ordinary impurity Anderson model and its known properties are recovered. For 
other positions, there is a partial destructive interference in the cotunneling processes involving the 
two doublets and as a consequence, the Kondo temperature and the height and width of the central 
peak (for bias voltage Vb near zero) of the differential conductance G = dl/dVb (where I is the 
current) are reduced. In addition, two peaks at finite Vb appear. We study the position of these 
peaks, the temperature dependence of G and the spectral densities. Our formalism can also be 
applied to carbon nanotube quantum dots with intervalley mixing. 

PACS numbers: 73.23.-b,73.22.-f, 75.20.Hr 



I. INTRODUCTION 

Single molecule electronic devices are being extensively 
studied because they offer perspectives for further minia- 
turization of electronic circuits with important potential 
applications 1-5 , and also due to its intrinsic interest in 
basic research, for example as realizations of the Kondo 
effect, which manifests experimentally by an increased 
conductance at low temperatures-' 7 . In addition, quan- 
tum phase transitions involving partially Kondo screened 
spin-1 molecular states were induced changing externally 
controlled parameters^—. These experiments could be 
explained semiquantitatively using extensions of the im- 
purity Anderson model treated with either the numer- 
ical renormalization group (NRG) 9 - - — or with the non- 
crossing approximation (NCA)i2r— . This approximation 
allows calculations out of equilibrium and in particular 
at finite bias voltage. 

The conductance through benzene- 1,4-dithiol has 
been measured using a mechanically controllable break 
junction 14 . More recently, other molecules containing 
benzene rings or related phenyl groups were studie d 15 ' 16 . 
The states nearer to the Fermi level in benzene are built 
from the 2p states of the C atoms which lie perpendicu- 
lar to the plane of the molecule (the so called 7r states). 
Transport through aromatic molecules and in particular 
through benzene in different geometries has been calcu- 
lated by several groups^ 7 - - — . Hettler et alJZ started from 
an exact calculation of an extended Hubbard model for 
the 7r states, and included the coupling to the leads per- 
turbatively (for couplings smaller than the temperature) , 
to calculate the current through neutral benzene under 
an applied bias voltage Vb for zero gate voltage V g , includ- 



ing radiative relaxation. Cardamone et alJ£, described 
the molecule by a Hubbard model supplemented by repul- 
sions at larger distance (the so called Pariser-Parr-Pople 
(PPP) Hamiltonian 22 ) using the self-consistent Hartree- 
Fock approximation. The conductance of the effective 
one-body problem for small V g has been calculated using 
the Landauer-Biittiker formalism—, The authors propose 
to exploit the destructive interference and to control the 
conductance of the device introducing symmetry break- 
ing perturbations. Similar results were obtained using 
ab initio calculations — , which have also the drawback of 
neglecting the effect of correlations. 

This effect together with that interference were in- 
cluded in more recent works which started from the ex- 
act solution of the PPP mode l 20 ' 21 , as well as in pre- 
vious works which discussed the conductance through 
molecules or rings of quantum dots (QDs) threaded by 
a magnetic flu»2i _ — . The conductance G depends on 
which sites of the molecule are coupled to the leads 
and on interference phenomena related to the symme- 
try of the syste m 20 ' 21 ' 29 . For certain conditions G can 
be totally suppressed and restored again by symmetry- 
breaking perturbations 21 . While the main conclusions 
are safe in general, the coupling to the leads was in- 
cluded in some perturbative approach, either by a Li- 
ouville equation method 20 , or using the Jagla-Balseiro 
formula (JBF) for the conductanc e 24 ' 26 , which assumes a 
non-degenerate ground state. These approaches miss the 
Kondo effect which is non-perturbative in the coupling 
to the leads 2 ^. This is particularly important when the 
gate voltage is such that the benzene molecule becomes 
charged, because the conductance is larg 

The JBF has also been used to predict that the trans- 
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mittance integrated over a finite energy windo w ' 25 i 27 
and the equilibrium conductanc e) 25 ' 29 through a ring of 
strongly correlated one-dimensional systems display dips 
as a function of the applied magnetic flux at fractional 
values of the flux quantum, due to destructive interfer- 
ence at crossings of levels with different charge and spin 
quantum numbers 2 ^ (as a consequence of spin-charge sep- 
aration). Recently, we have confirmed the presence of 
dips in the current under a finite applied bias voltage at 
low temperatures, using the low-energy effective Hamil- 
tonian, consisting of two doublets and a singlet, for the 
case of perfect destructive interference 3 ^. At equilibrium 
(zero bias voltage) , the conductance as a function of tem- 
perature for this particular case has been calculated with 
NRG 3 -. The interplay between interference and inter- 
actions were also studied for spinless electrons in mul- 
tilevel systems 3 ^, and benzene attached to two leads^ 3 -. 
Although at first sight they might seem artificial, spin- 
less models describe effective Hamiltonians for realistic 
systems under applied magnetic field a 34 ' 35 . 

Interference phenomena were observed in systems with 
QDs. In a multilevel QD, the crossing of levels and the 
ensuing destructive interference has been induced apply- 
ing magnetic field to a system with large, level-dependent 
g factors^ 4 -. Aharonov-Bohm oscillations were observed 
in systems involving two quantum dots (QDs)^ 6 .. 

In this work, we construct the low-energy effective 
Hamiltonian for a ring of n sites with one orbital per 
site and symmetry that includes the point group C nv (or 
C n ), weakly coupled to two conducting leads, retaining 
for n even, the lowest singlet with n particles and the two 
lowest doublets with n + 1 particles (electrons or holes 
depending on the sign of the applied gate voltage V g ). 
For n odd, the charge of doublets and singlets are in- 
terchanged. Using a gauge transformation, three of the 
four hopping matrix elements between the doublets and 
the leads can be made real. The phase of the fourth one 
<j) is in general different from zero and depends on the 
position of the leads and the wave vectors of the states 
involved. While the coupling to the leads should be small 
in such a way that the neglected states do not affect the 
calculations^ 7 -, we treat this coupling in a self-consistent, 
non perturbative way, using the the Kcldysh formalism 
of the NCA for problems out of equilibriu m 13 ' 38 ! 39 , ap- 
propriately extended for this problem. We calculate the 
non-equilibrium conductance for case of benzene (n = 6) 
in a regime of V g for which the doublets are favored, lead- 
ing to Kondo effect. 

For the case of one doublet, comparison of NCA with 
NRG results 40 , shows that the NCA describes accurately 
the Kondo physics. The leading behavior of the differ- 
ential conductance for small voltage and temperature 41 - 
agrees with alternative Fermi-liquid approaches 4 ^ 3 -, and 
the temperature dependence of the conductance practi- 
cally coincides with the NRG result over several decades 
of temperature 4 ^. A shortcoming of the NCA is that, at 
very low temperatures, it introduces an artificial spike at 
the Fermi energy in the spectral density when the ground 



state of the system without coupling to the leads (zero 
hybridization) is non-degenerate, although the thermo- 
dynamic properties continue to be well described 38 . This 
is not important in the present work, where the doublets 
are below the singlet and there is no applied magnetic 
field. Another limitation of the method is that it is re- 
stricted to temperatures above ~ T^-/20, where Tk is the 
Kondo temperature. In our case, this limitation is not 
important either because as we shall see, the conductance 
has already saturated at temperatures above this limit. A 
virtue of the NCA that becomes important in our case, 
is its ability to capture features at high energies, such 
as peaks in the spectral density out of the Fermi level, 
which might be lost in NRG calculations 44 . An exam- 
ple is the plateau at intermediate temperatures observed 
in transport through Cgo molecules for gate voltages for 
which triplet states are importan^^, which was missed 
in early NRG studies, but captured by the NCAHii 3 .. 
More recent NRG calculations using tricks to improve 
the resolution^ 5 ., have confirmed this plateau 1 - -. 

For the particular case of perfect destructive interfer- 
ence <f> = it (which does not correspond to benzene), the 
spectral densities 4 ^ and the non-equilibrium current^ ., 
were calculated before for symmetric coupling to the 
leads. In this case, the model is equivalent to an SU(4) 
impurity Anderson model, used by several authors 4 ^— 
to interpret transport experiments in C nanotube o 51 ' 52 , 
and more recently in Si fin- type field effect transistors 53 . 
If in addition a finite splitting 6 between the doublets is 
allowed, the symmetry is reduced to SU(2). The transi- 
tion between SU(4) and SU(2) has been studied by sev- 
eral authors using the NC A 30 i 46 ' 49 i 53 . In particular we 
found that the reduction of the Kondo temperature with 
S agrees very well with a simple formula obtained from a 
variational calculatio n 30 ! 46 . An important difference be- 
tween our model and those used for other systems is the 
connectivity with the leads, which imposes that in our 
case, the equilibrium conductance G(V& = 0) is zero in 
the SU(4) limit 6 = 0, and therefore, it is not propor- 
tional to the total spectral density of states^ -. For (j> = it 
and any 5, it is possible to calculate G(0) in terms of 
static occupations, using a generalized Friedel-Langreth 
sum rule^ -, due to the simplicity of the "orbital" field 
6, which acts in the same way as a magnetic field, and 
allows to relate the phase shift for each channel and spin 
with the respective occupancies. In turn, these phase 
shifts enter a general Fermi-liquid expression for the con- 
ductance obtained by Pustilnik and Glazman 54 . This is 
not possible for general <fi. Therefore a novel treatment is 
needed to calculate the conductance even at equilibrium. 

Recently, it has been shown that in carbon nanotubc 
QDs with disorder induced valley mixing, the SU(4) sym- 
metry is broken and the tunnel couplings to metallic leads 
become complex and depend on the applied magnetic 
field^ 5 -. The present formalism is also suitable to study 
this case. 

The paper is organized as follows. In Section [TT] we de- 
scribe the construction of the effective Hamiltonian. The 
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FIG. 1. (Color online) Scheme of the relevant matrix elements 
at low energies. 

approximations and the equation for the current are pre- 
sented in Section ITTT1 Section llVl contains the numerical 
results for transport through the benzene molecule and 
their interpretation. Section [V] contains a summary and 
a short discussion. In appendix 1X1 we show that the NCA 
conserves the current, extending the existing demonstra- 
tion for the case of one doublet^ 8 .. Appendix [B] contains 
some numerical details. 



II. MODEL 

The effective model, which is represented at the right 
of Fig. [TJ contains a singlet with total wave vector Kq 
(usually or w) and two doublets with wave vectors K\ 
and K 2 , that represent the states of lowest energy of two 
neighboring configurations of a ring of n sites with sym- 
metry C nv . In the case of benzene, they correspond to the 
singlet ground state, invariant under rotations ( Kq = 0) 
and two degenerate doublets with total wave vectors ±K, 
which are the ground state of the molecule for one added 
electron or hole, depending on the sign of the applied gate 
voltag o 20 ' 21 ' 29 . For one added hole K = tt/3, while for 
one added electron K = 2ir/3. The effective Hamiltonian 
is 

H = E„\s){s\ +'^2E i \i<r)(ia\ + e vk & vk(T c vk< , 

id uka 

ivka 

where the singlet \s) and the two doublets \ia) (i = 1,2; 
a =t or I) denote the localized states, c\ ka create con- 
duction states in the left (y — L) or right (y — R) lead, 
and V" describe the hopping elements between the leads 
and both doublets, assumed independent of k. This hop- 
ping element is calculated as 2 ^ 

V? = t u {ia\c) v<r \8) y (2) 

where c\j CT creates an electron (or hole) with spin a at the 
7T orbital at site j (1 to n) of the molecule, j v denotes 
the site connected to the lead v, and t v is the hopping 
between this site and the lead v. 

Choosing adequately the phases of the gauge transfor- 
mation \ia) — > e lVi \ia) : both can be made real and 



positive, and by reflection symmetry Vi — V 2 . Using 
rotational symmetry, it is easy to see that 2 ^ 



V t R = (t R /t L )V^ exp [-i(j R - j L )(Ki - K )\ , (3) 

where Ki is the wave vector of \ia). The phase of 
can be absorbed by a gauge transformation in the CRk a , 
rendering it real and positive. The remaining matrix el- 
ement is V 2 R = Vfe-**, where (j) = (j R - j L )(K 2 - K x ). 
This result is general for a ring geometry. 

For benzene K 2 — K\ = ±27r/3. Then, if the leads are 
connected in the para position {jr—Jl = 3) , <p = Q, while 
in the ortho {jR—jh = 1) or meta (jR—jL = 2) positions, 
4> = ±2tt/3. Note that the sign of <p does not affect our 
results. For simplicity we shall assume |t«/tL| = 1 in 
the calculations presented here. Then, the hopping of 
the leads to the relevant states of the benzene molecule 
are Vf = V 2 L = = V, V 2 R = Ve' 1 *, with = 
for the para position, and (f> = ±2ir/3 for the other two 
possibilities of connecting the leads. 

For = 0, the state \Ba) = - \2a))/V2 decou- 

ples from the leads and the transport properties of the 
system are the same as those for a single level [\Aa) = 
(|1(t) + |2cr))/V2] QD connected to he leads, with hop- 
ping v2V, which are well know n 38 1 56 " —. Another known 
limit of the model, which cannot be realized in benzene 
molecules, but in rings of a number of sites multiple of 
four— is cf) — it. In this case, \Aa) (\Ba)) is coupled 
only to the left (right) lead and the model is equivalent 
to an SU(4) Anderson mode l 30 ' 46 . However, in contrast 
to the case of C nanotubes^r— and Si fin-type field ef- 
fect transistors^ 3 -, the conductance vanishes identically 
as a consequence of perfect destructive interferenc e 30 ' 31 . 
An important technical difference is that in our model, 
the hybridization matrices of the states \ia) with the left 
and right leads are not proportional for <\> ^ 0. If the 
degeneracy between the levels can be lifted (for example 
applying an external flus 2 ^) , the symmetry is reduced to 
SU(2) and a finite conductance is restore d 30 ' 31 ' 46 . In C 
nanotubes with disorder induced valley mixing 5 -§., the ef- 
fective model becomes very similar to ours, including a 
finite level splitting and a phase (j>, which depends on the 
applied magnetic field. 

We shall use our previous results for the case c/> = n 
and a finite level splitting 6 = E 2 — E\ to help in the 
analysis of the results presented here for benzene. 



III. THE FORMALISM 

In this section we describe the extension of the non- 
crossing approximation (NCA) applied before for the 
Anderson model with infinite on-site repulsion out of 
equilibriu m 38 ' 39 , to our effective Hamiltonian. 
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A. Representation of the Hamiltonian with slave 
particles 

An auxiliary boson 6, and four auxiliary fermions fi a 
are introduced, so that the localized states are repre- 
sented as 



|s>=&t| ),|«T> = /t,|0>, 



(4) 



where |0) is the vacuum. These pseudoparticles should 
satisfy the constraint 



(5) 



propagators by the dressed ones, the expressions for the 
lesser self-energies take the form 

n<( W ) = - r r m / ^(i-^-w))^), 

ulma 

^»=-J2 r imj ^f V {u-u')D<{u>), (8) 



where 



r&( W ) = 27r£^%-<*fc) 



(9) 



Introducing it by a Lagrange multiplier, the effective 
Hamiltonian takes the form 



H' = (E s + \)tfb + + X)flf ia + J2 **4a<W 

icr vkcr 

+ J2(vfflbc k „ a +K.c). (6) 



The NCA solves a system of self-consistent equations to 
obtain the Green functions of the pseudoparticles (de- 
scribed below), which is equivalent to sum an infinite 
series of diagrams (all those without crossings) in the 
corresponding perturbation series in the last term of H', 
and afterwards project on the physical subspace of the 
constraint. The main new difficulties compared to the 
case of one doublet arise as a consequence of the matrix 
structure of the pseudofermion Green functions and self 
energies, which are absent in the SU(4) case or in this case 
with simple SU(2) symmetry breaking perturbations, as 
a magnetic^ 3 - or " orbital"^! d& field. 



B. Green functions 

The lesser and greater Keldysh Green functions for the 
psedoparticles for stationary non-equilibrium processes 
are defined a o 60 ' 61 



D < (t-t') = 
Gf^(t-t') = 
D>(t-t') = 



-<(6t(f)6(t)), 
-iifiaWlV)), 



(7) 



the retarded and advanced fermion Green functions are 
= miGlAV + G^Jt)}, G^ a = GV j<a +Gf ji(T - 
Gfj a , and similarly for the boson Green functions replac- 
ing Gij^ by D. 

These Green functions correspond to the interact- 
ing (dressed) propagators and are determined selfconsis- 
tently within the NCA. Evaluating the corresponding di- 
agrams in second order in the V" , and replacing the bare 



assumed independent of ui. Similarly, the greater self- 
energies become 

n>( w )= r imJ ^M^-^g^M)- (io) 

vlrarr 

As in the case of one doublet only^, in the expressions 
for the retarded self energies, E r (f) = -0(t)(£<(i) - 
E > (t)), E < can be neglected in comparison with E > due 
to the effect of the constraint. Then after Fourier trans- 
forming one obtains 



v.- , ■ ; duj> S ^>( w/ 



2tt u> — oj' + ir\ ' 

du)' n>(w') 



UJ — UJ 1 + IT) 



(11) 



where r\ is a positive infinitesimal. The advanced self 
energies E^ a and II a are obtained changing the sign of 
j] in the expressions above. 

For the case of SU(4) symmetry or when this symmetry 
is broken by simple symmetry breaking field o 30 ' 31 i 46 ~ 50 i 53 
((ft = it and symmetric coupling to the leads in our ef- 
fective model), the fermion Green functions and self 
energies are diagonal in an appropriate basis and the 
self-consistency problem simplifies considerably. In the 
general case (including <f> = ±2ir/3 for benzene) one 
has to solve a matrix Dyson equatio n 60 i 61 which in- 
cludes not only the Keldysh (or +, - branch index in the 
Keldysh contour), but also the doublet index i = 1,2 
in the fermion case. For the retarded fermion Green 
functions and self energies, combining G\a a in a 2 x 2 
matrix G r and similarly for self energies and unper- 
turbed Green functions CT = Sij(oj — Ei)^ 1 , with 
Ei = E — i + A the Dyson equations take the simple 
form G r = g r + g r E r G r . Solving the system for the 
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G\j a we obtain 





-(uj- E 2 
D 


^22 




U 12,<t — 


g( S 12,J, 






^21,(7 — 








U 22,a — 




^11 





(12) 



where 



El - Sii CT )(w — E 2 — 



-^12,(7^21,. 



For the boson, which has no subscripts, one has 



E. - IF 



(13) 



The advanced Green functions can be obtained from 
the retarded ones by the replacement rj — > —rj. 

The remaining equations that relate the lesser and 
greater pseudoparticle Green functions with the corre- 
sponding self energies are (in compact matrix form for 
the fermions)^ 2 - 



G? 



G r S^G Q , 

D r n^D a . 



(14) 

We solve numerically the system of integral equations 
([8]) to ([14]) for the pseudoparticle Green functions and 
self energies. The details are given in appendix IB1 

For the calculation of the current, one needs the Green 
functions of the physical fermions d\ a = \io)(s\ = fj a b. 
These functions, which we identify with the subscript d, 
are defined in the same way as the pseudofermion ones 
[see Eqs. (J7J], replacing f icr by d icr . In terms of the 
auxiliary-particle Green functions, the lesser and greater 
physical Green functions take the forrn^ 

W -G< ji<T ( U ' + U )D>( U '), 



2nQ 
duj' 
2^Q 



GZJu/ + u>)D<(tS), (15) 



where Q = (b^b + j\ a Ua) for a given Lagrange mul- 
tiplier A (see appendix |B|) . 



where the + (-) sign correspond to the case in which the 
doublets have one more electron (hole) than the singlet, 
f v (w) = [exp[(w— p^/kT]^!]^ 1 where p v is the chemical 
potential of the lead v, and T v , and are 2 x 
2 matrices with matrix elements given by Eqs (|9]) and 
(1151) . In particular, taking the unperturbed density of 
conduction states per spin p = 1/(2D) where 2D is the 
band width, and symmetric coupling to the leads, we 
have 



r L 



nV 2 ( 1 1 
D \ 1 1 J ' 

irV 2 t 1 
D { e-^ 1 



(17) 



Note that unless = 0, T L is not proportional to T R . As 
a consequence, the trick used to relate the conductance 
at Vb = with the spectral density of state o 64 i 65 cannot 
be used. We calculate the conductance G = dl/dVb by 
numerical differentiation of the current even for Vb — > 0. 
The traces appearing in Eq. (|16l) have the form 



Tr(T R G>; 



ttV 2 

IT 

+i 



[Gfi 



G g 
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cosW>)(G> 



Gf 2 ) 



sin(^)(G> 



Gf 2 ), 



(18) 



and similarly for Tr(r L Gg") replacing <j> by 0. 

Note that from the definition of the Green functions 
[see Eq. fl7J|], one realizes that the complex conjugates 
of the lesser and greater Green functions satisfy G?* (t) = 
-G<l(—t) and after Fourier transform Gfj(ui) = -G^(w). 
Then, G|(w) and G>(cj) + G>(w) are pure imaginary, 
while Gf^u) — Gf 2 (i-u) is real. 

In appendix A we show that the current from the left 
lead to the molecule equals that from the molecule to the 
right lead, so that the current is conserved within the 
NCA. Other approximations, like perturbation theory in 
the Coulomb repulsion out of the electron- hole symmetric 
point do not conserve the current^, unless some tricks 
are used§2i£L£S. 



IV. NUMERICAL RESULTS 



C. Equation for the current 

Using general formulas for the current through a region 
with interacting electron o 64 i 65 and the relation G r — G a = 
G > — G < between Green functions, the current of our ef- 
fective model for a spin degenerate system (without ap- 
plied magnetic field) can can be written as 

i = ±jj ^Tr[(r L /LH-r R /flH)G>H 
+ ((1 - / L M)r L - (1 - /*H)r R )G<M], (16) 



For the numerical calculations, we assume a constant 
density of states per spin of the leads p between — D 
and D. We take the unit of energy as the total level 
width of both doublets: T = Tf t + r£ = 4irpV 2 , i = 
1,2. Also r = 2A, where A, called the resonance level 
width, is half the width at half maximum of the spectral 
density of states of particles for each level in the non- 
interacting case. We restrict our present study to gate 
voltages such that E s — Ei 3> A, for which correlations 
play a more important role, the Kondo effect develops 
and the conductance is higher. Without loss of generality, 
we take 6f = E 8 = 0, where cf is the Fermi level of 
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the leads without applied bias voltage. We define Ed = 
Ei, 8 — E 2 — Ei. For a bias voltage the chemical 
potentials are (1l = eVb/2, = — eVb/2. 

In the numerical results an important low-energy scale 
is the Kondo temperature Tk- It is known that for im- 
purity Anderson models like ours for 8 = in the Kondo 
regime, the spectral density shows one charge-transfer 
peak at energy Ed below the Fermi energy of total width 
of the order of AT for SU(N) models^^. For finite on- 
site Coulomb repulsion U there is another charge-transfer 
peak at energy Ed + U, but this is shifted to infinite ener- 
gies in our case. In addition, there is a narrow peak near 
the Fermi energy of width of the order of Tk- In this 
work we define Tk as the half width at half maximum of 
this peak. For finite 8 additional peaks appear as shown 
in Section HVB 



A. Conductance out of equilibrium 

In Fig. [5] we show the differential conductance G = 
dl/dVb as a function of bias voltage V&, for the leads 
connected at 180 degrees (in the para position), for which 
(j) = 0, and in the other two positions (<j> = ±27r/3), for 
several values of the energy of the doublets Ei = E 2 = 
Ed- The charge-transfer energy tp — Ed is the energy 
necessary to take a particle (electron or hole) from the 
Fermi energy, and add it to the localized singlet, forming 
a doublet, in absence of the hopping to the leads. It is 
tuned modifying the gate voltage V g : Ed decreases with 
increasing V g . The temperature in the curve was fixed at 
T = 0.057V. Since we have assumed symmetric coupling 
to both leads, G(—Vb) — G(V&) and only positive 14 are 
duisplayed in the figure. 

For the para position, as explained in Section [2] the 
problem is equivalent to the transport through a single 
level quantum dot, and G shows a single peak centered 
near zero voltage and width of the order of 2fcTx/e 69 . 
We remind the reader that half the width of G(H) for 
T = times the electric charge e, half the width of the 
spectral density and the temperature (times the Boltz- 
man constant) T at which G(T) for T4 = falls to half of 
the maximum value, are all of the same order—. This is 
natural from the expected universal behavior of the (sin- 
gle level) Anderson model in the Kondo regime, in which 
only one energy scale Tk survives. This Kondo tempera- 
ture varies exponentially with the charge-transfer energy 
Ed, and as a consequence, the width of the central peak 
observed in the figure also decreases exponentially with 
increasing \Ed\, as observed in Fig. [2] 

When the molecule is attached to the leads at 60 or 
120 degrees, three main changes occur in G(Vb) in com- 
parison to the previous case: i) The maximum of G(0) 
is lower, ii) The width of the central peak is narrower 
and scales in a different way with Ed, and iii) two new 
peaks appear at finite V&. i) is due to the effect of partial 
destructive interference between the current transmitted 
by the two doublets. In the para position, for which all 
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FIG. 2. (Color online) Differential conductance at low tem- 
perature through a benzene molecule as a function of bias 
voltage for several values of the energy of the doublets, with 
the leads connected in the para (red dashed line) and ortho 
or meta (full black line) positions. 



hybridization to the leads have the same sign (cf> = 0), 
G(0) equals the quantum of conductance 2e 2 /h for sym- 
metric coupling to the leads (as we have assumed in the 
present calculations). In the opposite case (hypothetical 
for benzene) of 4> = ir, there is perfect destructive inter- 
ference and G vanishe d 30 ' 31 . Transport through the ben- 
zene molecule connected in the ortho or meta positions, 
correspond to an intermediate situation (<j) = ±27r/3), 
and therefore, an intermediate value of G(0) is expected. 
For asymmetric coupling to the leads, as it is usually the 
case in transport through molecules in devices built by 
electromigration&i£, G(0) decreases strongly and it is not 
possible to distinguish between different positions of the 
leads connected to the benzene molecule from the maxi- 
mum value of G. Note that this value decreases slightly 



7 




0.08 



FIG. 3. (Color online) Differential conductance for <j> — n and 
different values of the level splitting. Other parameters are 
E d = -4,T = 0.05T K . 



with decreasing Tk- The feature, which is probably the 
most useful to distinguish experimentally among the dif- 
ferent positions of the conducting leads is the presence 
of the side peaks in G(Vb), which seem to depend weakly 
on the charge transfer energy tp — Ed, in contrast to the 
width of the central peak. 

In order to help in the interpretation of the results, 
we calculate G(Vb) for <p = it, but allowing for a finite 
splitting 5 = E% — E\ between both doublets. S acts as a 
symmetry breaking field on the SU(4) Kondo effec t 30 ' 46 . 
The result is shown in Fig. [3j Clearly, G(Vb) is quali- 
tatively similar to the corresponding result for benzene 
with leads connected in the ortho or meta positions. The 
conductance, which vanishes in the SU(4) limit S = is 
restored by a finite S and two peaks at eVb = ±<5 appear. 
This similarity suggest to interpret the results for ben- 
zene starting from those for 4> = n, S = and thinking the 
difference between the coupling for = ±27r/3 and 
4> = 7T as a perturbation. This perturbation, in second 
order, introduces among other effects, an effective mixing 
between the levels, proportional to V 2 /\Ed\, which leads 
to a splitting of the doublets. In fact, the position of the 
satellite peaks in Fig. [2] (ranging from 0.3 for Ed = —3 
to 0.22 for Ed = —5) is roughly consistent with a l/|£d| 
dependence. 

Note that in the case cj) — in which this perturba- 
tion is expected to be the largest, only the symmetric 
combinations of left and right lead states hybridize with 
the impurity, while the antisymmetric ones remain decou- 
pled. The spectral density of the latter is a delta function 
6(uj — Ed), without weight near the Fermi energy. This 
makes clear that actually the effective hopping of the two 
resulting split doublets is different and the proposed in- 
terpretation is rather qualitative. 




Ed: 



-0.5 -0.4 -0.3 -0.2 -0.1 Oq) 0.1 0.2 0.3 0.4 0.5 

FIG. 4. (Color online) Equilibrium spectral density for the 
benzene doublets as a function of frequency in the ortho and 
meta positions for different temperatures. 




FIG. 5. (Color online) Total spectral density at equilibrium 
as a function of frequency for the effective model with <j> = n, 
T = 0.05Ta' and different values of the level splitting. 



B. Spectral densities 

As a further test of the interpretation outlined above, 
we compare the spectral density of states of both models. 
The spectral density 

p icr (L>) = (G dtia (uj + irj) - G d ua{u - irj))/(2iri) (19) 

of both doublets for benzene connected to the leads in 
ortho or meta positions is represented in Fig. 2] for dif- 
ferent temperatures and energies near the Fermi energy 
(the charge-transfer peak 4 ^ is not shown). As T decreases 
below Tk, not only a peak develops near the Fermi en- 
ergy, but also two side peaks (the most prominent for 
positive frequencies) are clearly present. In Fig. [S] we 
show the total spectral density for <p = tt, very low tem- 
peratures and different values of S. In this case, the peak 
near to the Fermi energy corresponds to the lowest dou- 
blet and the peak at energy near S corresponds to the 
highest doublet^,. In contrast, each peak in Fig. 2] is 
expected to come from some linear combination of both 
doublets which is not easy to identify. Another difference 
apparent in the figure is that the side peak at positive 
frequencies for benzene is sharper and more asymmetric 
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FIG. 6. (Color online) Conductance as a function of temper- 
ature for leads connected in the ortho or meta position and 
several values of the energy of the doublets. 

than the corresponding one for <fi = it and finite 5. This 
might be due in part to a smaller effective hybridization 
of the excited doublet in the case of benzene, leading to 
a narrower peak, as in the limit <f> — > discussed above. 

In spite of these differences, the spectral density for 
cj> = 7T and finite 5 shows the development of a satellite 
peak departing from the central one as S increases. The 
spectral density for benzene shows the same qualitative 
features (a peak near the Fermi energy and another one 
at positive frequencies) which can be interpreted, in a 
qualitative first approximation, as coming from an effec- 
tive S. This peak at finite energy in turn gives rise to the 
side peaks in G(Vb) (although the spectral densities are 
modified with the applied bias voltage and the peaks are 
blurred). 



C. Temperature dependence of the conductance 

In Fig. |51 we show the evolution with temperature of 
the conductance G(T) for Vb = 0, several values of Ed, 
and conducting leads in the ortho or meta position. G(T) 
was obtained differentiating the current, as explained in 
Section [ill CI In order to be able to display all curves in 
the same scale, a logarithmic scale of relative tempera- 
tures T/Tk is used, where the characteristic energy scale 
Tk is given by the half width at half maximum of the 
equilibrium spectral density pi a {uj). The resulting value 
of T K is 3 x 10" 2 , 3.6 x 1CT 3 and 4.8 x 1(T 4 for E d = -3, 
-4 and -5 respectively. The peak at finite energy in the 
spectral density at low temperatures (which lies between 
0.3 for Ed = —3 to 0.22 for Ed = —5) results in a struc- 
ture in G{T) at T/T K near 10, 70 and 400 for E d = -3, 
-4 and -5 respectively. While for Ed = — 3 this structure 
is hidden in the main peak, a kind of plateau is evident 
for the other values of Ed at the corresponding positions, 
as shown in the inset of Fig. [5] 

This structure is reminiscent of the plateau observed 
by transport experiments in the "triplet" side of the 
quantum phase transition in Cgcr^^, and explained 



by NC A 10 ' 12 i 13 and NRG calculations with improved 
resolution^ ., In our case, to represent accurately the high 
energy features, it was necessary to use a special mesh 
in the frequency axis with more points near the high en- 
ergy peak in the pseudofermion spectral densities (see 
appendix [B]) . 

In previous wor k 30 ' , we have found that for <fi = it 
and any level splitting 6 = E2 — E\, Tk (except for a 
constant of order unity) is very accurately given by the 
expression 

T K = {{D + 5)Dexp[TrE 1 /(2A)}+6 2 /4} 1/2 -6/2, (20) 

obtained from a simple variational wave function that 
generalizes the proposal of Varma and Yafet21 for the 
simplest impurity Anderson model. This equation inter- 
polates between the SU(4) limit S = 0, for which Tk has 
the largest value T K = Dexp [7r£'i/(4A)] and the one- 
level SU(2) limit 5 -> +00. The effect of 6 on T K is small 
when S < T K , while for larger 5, Tk decreases strongly. 
Note that T K coincides with the Kondo temperature of 
the equivalent one-level effective Anderson for benzene 
connected in the para position (because the effective one- 
level resonant level width contains a factor 2, see Section 
ITT|) . Therefore, we can interpret the fact that the effective 
energy scale in the ortho and meta positions is smaller 
than that in the para positions (compare the widths of 
the central peaks in Fig. [2]) as an effect of the effective 
level splitting caused by a phase <f> different from it. 

In other words, it seems that the position of the side 
peaks in G(Vb) (or spectral densities) and the plateau in 
G(T) indicate an effective S, which decreases the charac- 
teristic energy scale Tk ■ While we do not expect Eq. (|2"0")) 
to describe accurately Tk using this effective S, the fact 
that the difference between the energy scales for the para 
and the other positions is more noticeable for S larger 
than T K {Tk for the para case), suggests that this equa- 
tion might be useful for a qualitative understanding, and 
supports the interpretation of our results. 

V. SUMMARY AND DISCUSSION 

We have constructed the effective Hamiltonian for 
transport through a symmetric ring (point group includ- 
ing C nv as a subgroup) with one orbital per site, includ- 
ing a singlet and two degenerate doublets of a neigh- 
boring configuration for the isolated ring. This includes 
for example the ground state of the isolated ring for one 
electron per site, and another with one added electron or 
hole. The resulting effective generalized Anderson Hamil- 
tonian describes however more general situations. An 
extension to the case in which a magnetic flux is added 
(reducing the symmetry to C n and breaking the doublet 
degeneracy) is trivial. Partial destructive interference oc- 
curs in general when two levels with N particles are near 
to another one with iV±l particles^ 9 -. Therefore, the situ- 
ation of one singlet and two doublets appears frequently. 
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Prom the wave vectors of these states, it is easy to calcu- 
late the phase <f> of the effective model following the lines 
of Section [TTJ Therefore, we expect that our formalism 
can be used in a variety of physically relevant systems. 

We have used the resulting effective Hamiltonian to 
calculate the non-equilibrium transport through a ben- 
zene molecule with conducting leads connected in differ- 
ent positions, using an appropriate generalization of the 
NCA to include the Kondo effect as well as effects of par- 
tial destructive interference between the transport chan- 
nels through both doublets. While for leads connected 
in the para position, the conductance G as a function of 
voltage and temperature is similar to that well known for 
the case of only one doublet, in the other positions the 
peak in G(Vb) near zero bias voltage \% is narrower and 
of smaller intensity as an indication of the relative posi- 
tion of the conducting leads. In addition, G(Vb) displays 
two additional side peaks at finite ±Vb, and a character- 
istic plateau is present in the conductance as a function 
of temperature G(T). These finite energy features are 
probably easier to detect experimentally as an indication 
of the relative position of the conducting leads. 

These results for the leads connected in the ortho and 
meta positions, which are due to partial destructive in- 
terference, can be interpreted starting from those with 
total destructive interference (corresponding to a phase 
4> = 7r) for which the effective model has SU(4) symme- 
try, and adding a symmetry breaking splitting S between 
effective doublets caused by the remaining term treated 
as a perturbation. 

For our results to be valid, the coupling to the leads 
should be sufficiently small, so that excited states, not 
included in the effective Hamiltonian, do not play an im- 
portant role. The effect of these states can be estimated 
for each particular case. See for example Ref. 37 . The 
magnitude of the coupling can be controlled for example 
in break junctions&2^. 

Our effective model, including a finite level split- 
ting S from the beginning, and a general phase 4> de- 
scribes carbon nanotube QDs with disorder induced val- 
ley mixing^. Therefore, our formalism supplemented by 
realistic calculations of 8 and <j> might be used to calculate 
the conductance of these systems. 

In this work, we have neglected the effect of phonons, 
which can modify the effective parameters of the 
mode l 72 i 73 , and also affect the interference phenomena 74 . 
We have also neglected electrostatic interactions between 
the leads and the molecules. In the weak coupling regime, 
it has been shown that interaction with image charges 
for the doped molecules leads to a symmetry breaking 
and a splitting of the degenerate levels that also affects 
the interference phenomena 7 ^. As we have shown, in the 
specific case of benzene, the coupling to the leads already 
originates a splitting of this kind. Therefore, we expect 
that these electrostatic effects are more relevant in the 
case of perfect destructive interference. 
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Appendix A: Current conservation 

Here we present a brief demonstration that within 
NCA, the current established between the left metallic 
contact and the central region (ring) 7j,, is the same as 
the current flowing from this region to the right metallic 
contact Ir. For simplicity we assume that there is no 
applied magnetic field. The demonstration can be easily 
extended to the case of inequivalent spins. Adding both 
spins the expressions of the currents are^ 

2ie f 

h = ±— / clw Tr{T L [/ L HG>H 



+ (i-/ L H)G<H]} 

-He 



(Al) 



Ir = ± 



, J duTv{T R [f R (io)G>(u) 
+ (I-/ i? H)G<H]}, (A2) 

where the symbol ± and the meaning of the matrices is 
explained in Section Mil [actually Eq. (fl~6| is the average 
of Eqs. ((SU and (jA2|)]. 

We use Eqs. (fT5|) to replace the physical Green func- 
tions by their expressions in terms of the Green functions 
for the auxiliary particles. Denoting G the 2x2 matrix 
with the pseudo fermion Green functions for a given spin,, 
one obtains 



-f L (u)D < (uj') Tr T L G (uj' + uj 



(A3) 



In = ±— / did 



dJ f R {u)D < {J) Tr \t R G > {u/ + u) 



h J J 2ttQ 
+ (l-f R {u))D>(u') U:[r R G < (u , +w) 



(A4) 



Using Eqs. dTU|) for the boson self-energies, and some 
algebra, the difference between both expressions becomes 

JW* = x J ^f[D < (u')n > (u>')-D>(u,')n<(u>>)], 

(A5) 

which is easily seen to vanish using D> = D r H>D a [see 
Eqs. CSM- 



Appendix B: Numerical procedure 

In this appendix we give some details about the nu- 
merical solution of the system of integral equations. 
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We first solve the system that determines the retarded 
and greater Green functions and self energies for the aux- 
iliary particles, Eqs. (fTU|) to (Q~4|). After the system con- 
verges, we solve self-consistenlty Eqs. © and (|T4"| for 
the lesser quantities. 

Due to the gauge symmetry of the Hamiltonian in 
the representation of the auxiliary particles, fi a — > 
exp(iAoi)/icr, b — > /j CT exp(iAo£)&, the Lagrange multiplier 
A in Eq. (IT5l) can be shifted by an arbitrary number Ag 39 -. 
This is equivalent to a shift in the energies of the slave 
particles or of the frequency as uj — > uj + Ao . Therefore, 

we can either set Q = (b^b + J2i<j flafiv) = 1 an d calcu- 
late A at each iteration or fix A, calculate Q and divide 



by Q in the evaluation of physical properties 38,39 . 

In the present problem, the spectral functions of the 
pseudo-fermions have two different peaks in contrast to 
the single one present in the ordinary Anderson model. 
We found that fixing A as the solution of the equation 
KeD(u, A) = at each iteration, the position of the low- 
est peak of the auxiliary fermion spectral densities coin- 
cide with the zero of energy (the Fermi level) With this 
chice of A, in order to make the numerical evaluations 
efficiently, we use a non equidistant frequency mesh 
with two different regions logarithmically spaced. One 
of them is aimed at the zero of energy uj — and the 
other one is centered near to the position of the second 
peak of the fermion spectral functions, u) = 8. Finally, to 
generate the non equidistant mesh we follow reference 3 - 9 -. 
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